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ABSTRACT 



Adult male fathead minnows (Pimephales promelas) were exposed to effluent from the City of Boulder, 
Colorado wastewater treatment plant (WWTP) under controlled conditions in the field to determine if 
the effluent induced reproductive disruption in fish. Gonadal intersex and other evidence of reproduc- 
tive disruption were previously identified in white suckers (Catostomus commersoni) in Boulder Creek 
downstream from this WWTP effluent outfall. Fish were exposed within a mobile flow-through expo- 
sure laboratory in July 2005 and August 2006 to WWTP effluent (EFF), Boulder Creek water (REF), or 
mixtures of EFF and REF for up to 28 days. Primary (sperm abundance) and secondary (nuptial tuber- 
cles and dorsal fat pads) sex characteristics were demasculinized within 14 days of exposure to 50% and 
100% EFF. Vitellogenin was maximally elevated in both 50% and 100% EFF treatments within 7 days and 
significantly elevated by 25% EFF within 14 days. The steroidal estrogens 17p-estradiol, estrone, estriol, 
and 1 7a-ethynylestradiol, as well as estrogenic alkylphenols and bisphenol A were identified within the 
EFF treatments and not in the REF treatment. These results support the hypothesis that the reproductive 
disruption observed in this watershed is due to endocrine-active chemicals in the WWTP effluent. 

© 2011 Elsevier BV. All rights reserved. 



1. Introduction 

Exogenous endocrine-active chemicals (EACs) can interact 
directly with endocrine regulatory systems (Tabb and Blumberg, 
2006), leading to the disruption of reproductive development and 
function in exposed organisms (McLachlan, 2001). Anthropogenic 
EACs enter aquatic environments through diverse human sources 
including military (Bordeleau et al., 2008), medical (Thornton et al., 
1996), industrial (Parrot et al, 2006), agricultural (Horrigan et al, 
2002), energy (Lye, 2000), and wastewater (Aerni et al, 2004; 
Barber et al, 2006; Johnson et al, 2005; Schultz et al, 2010; 
Vajda et al, 2008) activities. The resulting presence of EACs in 
the environment can elicit adverse effects in exposed organisms, 
including plant-microbe symbioses (Fox et al, 2001 ), invertebrates 
(Oberdorstor et al, 2001), and vertebrates (Mills and Chichester, 
2005). 

Effluent discharges from wastewater treatment plants ( WWTPs) 
contain a complex mixture of synthetic and biogenic EACs (Aerni 
et al, 2004; Barber et al, 2006; Johnson et al., 2005). Many 
surface waters that receive WWTP effluent have detectable lev- 
els of steroidal and non-steroidal estrogens (Kolpin et al, 2002; 
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Johnson et al, 2005), androgens (Kolodziej et al, 2003) and neu- 
roactive pharmaceuticals (Schultz et al, 2010). In some cases, as 
with 17p-estradiol, exogenous EACs are structurally and function- 
ally identical to endogenous steroid hormones. Their addition to 
endogenous hormone loads can destabilize vertebrate endocrine 
signaling (Vajda and Norris, 2006). 

Reproductive disruption has been observed in fishes inhabit- 
ing estrogen-contaminated waters (Jobling et al, 1998; Mills and 
Chichester, 2005). Recently, gonadal intersex, reduced sperm abun- 
dance, and elevated plasma vitellogenin levels in male fish were 
identified in white sucker (Catostomus commersoni) in a WWTP 
effluent-dominated reach of Boulder Creek, Colorado (Vajda et al., 
2008; Woodling et al, 2006). Effluent discharged from the Boulder 
WWTP comprises up to 80% of streamflow during low-flow condi- 
tions, and rarely contributes less than 35% (Murphy et al., 2003). As 
a result, native fish are exposed to a complex mixture of EACs at all 
stages of their life cycle (Barber et al, 2006; Vajda et al, 2008). 

The objective of this study was to determine whether the 
reproductive disruption previously observed in wild fish could 
be induced by the Boulder WWTP effluent. Likewise, the lack of 
reproductive disruption in the upstream water was evaluated. The 
biological results were then compared to the predicted relative 
estrogenicity of the two waters based on measured EAC concentra- 
tions. This study improves upon previous work that involves in situ 
exposure of caged fish by controlling for parameters known to effect 
endocrine status (i.e. temperature, oxygen, biotic environment) 
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(Aerni et al, 2004; Giesy et al, 2003; Harries et al, 1997; Purdom et 
al, 1994). Conversely, previously published controlled studies that 
expose fish to transported effluent that may not accurately reflect 
in-stream EAC concentrations due to chemical transformation after 
sample collection (Robinson et al, 2003). In the present study, con- 
trolled on-site, continuous-flow fish exposure experiments were 
conducted in a bioassay laboratory (Barber et al., 2007; Harries 
et al., 1999; Rodgers-Gray et al., 2001). Because EACs and other 
organic contaminants are known to leach from plastic materials 
(Soto et al., 1991), the system was constructed using only glass, 
stainless steel, and Teflon® . We evaluate primary and secondary sex 
characteristics of adult male fathead minnows (Pimephales prome- 
las) exposed to WWTP effluent and reference water from upstream 
of the WWTP outfall under controlled conditions in the field. Con- 
centrations of EACs were measured in the WWTP effluent and 
upstream reference water to evaluate their occurrence and con- 
centrations relative to observed biological effects. 

2. Materials and methods 

2A. Study site 

The WWTP for the City of Boulder, Colorado was selected 
because reproductive disruption was observed among free-living 
white suckers in Boulder Creek below the WWTP outfall (Vajda 
et al., 2008; Woodling et al., 2006) and the hydrogeochemistry 
of Boulder Creek (including estrogenic EACs) is well characterized 
(Barber et al., 2006; Murphy et al., 2003). This WWTP has an aver- 
age discharge of 0.74 m 3 s _1 (17 million gallons a day) and can 
contribute from <40% of stream flow during spring runoff con- 
ditions to >75% of the flow during base-flow conditions in late 
summer. Sewage from the City of Boulder was treated using a 
combined trickling filter/activated sludge process with nitrifica- 
tion/denitrification and chlorination/dechlorination, and hydraulic 
retention time of ~12h. The mean annual concentrations (Jan- 
uary 2005-December 2007) of ammonia, nitrate, biological oxygen 
demand, and total suspended solids in the WWTP effluent were 
6.6, 11.8, 14.5, and 6mgL _1 , respectively (City of Boulder). Two 
source waters were used for the fish exposure experiments, WWTP 
effluent (EFF), and Boulder Creek water from approximately 500 m 
upstream from the WWTP outfall (REF). 

2.2. On-site fish bioassay laboratory 

Fish exposure experiments were conducted during July 2005 
and August 2006 in a mobile laboratory deployed at the Boulder 
WWTP. During all exposure experiments, temperature (22 ± 1 °C), 
lighting, diet, aeration (>85% saturation), and flow (200mLmin _1 , 
six daily volume replacements) were controlled. Contaminant con- 
centrations in WWTP effluents and surface water can fluctuate 
on various time scales, and this design reflects these fluctuations 
rather than maintaining a constant exposure concentration. All sur- 
faces in contact with test solutions were made of glass, stainless 
steel, or Teflon®. Water from the EFF and REF sites was continu- 
ously pumped through Teflon® tubing to the lab using stainless 
steel pumps (Grundfos SQE; Geotech Environmental Equipment, 
Denver, CO) into separate 200-L stainless steel holding tanks posi- 
tioned above the laboratory. The two water sources were thermally 
equilibrated, then flowed by gravity to stainless steel splitter 
tanks that distributed the water to 10-L glass aquaria housing the 
fish. 

23. Experimental design 

Reproductively stimulated 12-month-old adult male fathead 
minnows were obtained from the U.S. Environmental Protection 



Agency (Cincinnati, OH) for experiments in 2005 and from Aquatic 
Biosystems (Ft. Collins, CO) for experiments in 2006. For both exper- 
iments, 10 randomly collected fish were processed upon arrival as 
initial controls. The remaining fish were weighed, measured, and 
randomly distributed to exposure treatments. In 2005, the fish were 
exposed to 100% EFF, 50:50 EFF: REF, and 100% REF treatments with 
four exposure durations (4, 7, 14, and 28 days) yielding 12 treat- 
ment groups. In 2006, the fish were exposed to 100% EFF, 50:50 
EFF:REF, 25:75 EFF:REF, 10:90 EFF:REF, 5:95 EFF:REF, and 100% REF 
treatments with two exposure durations (14 and 28 days) yielding 
12 treatment groups. In both experiments, 5-6 fish were placed 
into each of the 1 0-L aquarium. Mean body weight and total length 
[mean ± standard deviation (SD)] of the fish at the start of each 
study did not differ between treatment groups. 

Fish were fed frozen brine shrimp twice daily at 08:00 and 
20:00 h. The feeding rate was 2% of body weight per day, based on 
mean body weight on day-0 and adjusted for any deaths and fish 
sampled during the studies. Surplus food and feces were removed 
daily. Animal care and handling was in accordance with the Insti- 
tutional Animal Care and Use Committee of the University of 
Colorado. 

2.4. Tissue sampling and analysis 

In 2005, two tanks of fish (n = 6-10) were sampled from each 
treatment at 4, 7, 14, and 28 days. In 2006, one tank offish (n = 5-6) 
was sampled from each treatment at day 14 and day 28. Fish wet 
weight (0.01 g) and total length (mm) were recorded and used 
to calculate the condition factor [(body weight (g)/total lengths 3 
(mm)) x 100]. Nuptial tubercle prominence (Smith, 1974; Smith, 
1978) was scored as: 1 = tubercles not visible, 2 = tubercles visible 
as white discs, 3 = tubercles prominent, and 4 = tubercles prominent 
and protruding sharply. Dorsal fat pad prominence was assessed 
as 1 = not visible, 2 = soft discolored tissue, 3 = spongy thickened 
tissue, and 4 = dorsal hump with spongy tissue. Blood was col- 
lected into heparinized capillary tubes from the caudal vein, kept 
on ice, and centrifuged for 5min at 3000 xg within 3h. Hemat- 
ocrit was recorded and aliquots of plasma were frozen and stored 
at -40 °C until analyzed for vitellogenin by homologous enzyme- 
linked immunosorbent assay (ELISA) using an anti-fathead minnow 
kit (Biosense, Bergen, Norway). 

Gonads, liver, and spleen were dissected and weighed (0.001 g), 
and the gonadosomatic index (GSI) [(gonad weight (g)/body weight 
(g))xl00], hepatosomatic index (HSI) [(liver weight (g)/body 
weight (g)) x 100], and splenosomatic index (SSI) [(spleen weight 
(g)/body weight (g))xl00] were calculated. Freshly dissected 
gonads were preserved in 10% neutral-buffered formalin until pro- 
cessed for histology (Presnell and Schreibman, 1 997), at which time 
a whole testis from each fish was prepared. At least 1 whole cross- 
sections from each testis were evaluated on a 1-5 scale by light 
microscopy for the amount of mature sperm present (Pawlowski 
et al., 2004a,b): 1= sperm absent, 2 = weak, sperm prominent in 
<25% of tubules, 3 = moderate, sperm prominent in 25-50% of 
tubules, 4 = strong, sperm prominent in 50-75% of tubules, and 
5 = very strong, sperm prominent in >75% of tubules. 

2.5. Water sampling and chemical analysis 

Weekly water grab samples were collected starting at day- 
from the EFF and REF inflows into the on-site laboratory. The 
EFF and REF samples were analyzed for EACs by solid phase 
extraction (SPE) and continuous liquid/liquid extraction (CLLE), fol- 
lowed by gas chromatography/mass spectrometry (GC/MS) and gas 
chromatography/tandem mass spectrometry (GC/MS/MS) analysis 
(Barber et al., 2000). Surrogate standards were added to unfiltered 
water samples prior to extraction to evaluate method perfor- 



Table 1 

Average concentrations (±SD) of endocrine active chemicals detected in the Boulder wastewater treatment plant ( VWVTP) effluent (EFF) and Boulder Creek upstream from the WWTP outfall (REF). Also shown are chemical specific 
in vitro and in vivo 17|3-estradiol equivalency factors used to calculate 17|3-estradiol equivalency quotients. 



Compound 



In vitro estradiol 
equivalency factor 3 



In vivo estradiol 
equivalency factor 1 



Source 5 



2005 REF, 
n = 6 

(ngL- 1 ) 



SD 



2005 EFF, 
n = l 
(ngL- 1 ) 



SD 



2006 REF, 
n = 6 
(ngL- 1 ) 



SD 



2006 EFF, 
n = 6 

(ngL" 1 ) 



SD 



Steroids 

17(3-Estradiol 

1 7a-Ethynylestradiol 

Estrone 

Estriol 

17a-Estradiol 

Other estrogens 

4-Nonylphenol 

Bisphenol A 

1 ,2-Dichlorobenzene 

1 ,4-Dichlorobenzene 

4-rerr-Butylphenol 

4-tert-Pentylphenol 

4-Nonylphenolmonoethoxylate 

4-Nonylphenoldiethoxylate 

4-Nonylphenoltriethoxylate 

4-Nonylphenoltetraethoxylate 

4-Nonylphenolmonoethoxycarboxylate 

4-Nonylphenoldiethoxycarboxylate 

4-Nonylphenoltriethoxycarboxylate 

4-Nonylphenoltetraethoxycarboxylate 

4-t-Octylphenol 

4-t-Octylphenolmonoethoxylate 

4-r-Octylphenoldiethoxylate 

4-r-Octylphenoltriethoxylate 

4-t-Octylphenoltetraethoxylate 

Estradiol equivalency quotients 
Estradiol equivalency quotient 



1 
1.2 

0.2 



5.3E-5 



1 

25 
0.4 



0.003 



In vitro 
In vivo 



1-3 


<0.2 




2.1 


1.0 


<0.2 




3.2 


3.2 




2,3 


<0.2 




1.2 


0.7 


<0.2 




<0.2 






2,3 


<0.2 




75 


28 


<0.2 




60 


75 


> 




<0.2 




3.4 


3.2 


<0.2 




<0.2 




|S 




<0.2 




0.1 


0.3 


<0.2 




<0.2 






3,4 


210 


230 


2000 


1300 


73 


34 


240 


100 


Q 
TO 

a 




33 


46 


42 


44 


18 


27 


10 


5.2 


c 




<10 




<10 




<10 




120 


79 


£2 




81 


140 


2100 


3000 


<10 




1200 


360 


c 

Q 

cr. 

n 

>< 




<10 




68 


73 


<10 




47 


43 




<10 




65 


82 


<10 




1.9 


2.2 




480 


1200 


9800 


5800 


68 


71 


1500 


770 


o 




150 


400 


5000 


4000 


89 


140 


440 


380 


1 




62 


160 


2000 


1800 


<50 




280 


240 




<50 




450 


730 


<50 




39 


31 


8 




270 


330 


31600 


9400 


370 


220 


50000 


6800 


Kj 




300 


510 


42000 


7300 


200 


110 


34000 


13000 


2 




<50 




1500 


690 


<50 




510 


240 


Nj 




<50 




890 


740 


<50 




210 


72 


Uj 




17 


29 


200 


120 


<10 




30 


14 


Nj 




17 


42 


290 


260 


<10 




22 


11 


>-* 




82 


220 


2100 


2700 


8.7 


11 


53 


67 






<10 




260 


350 


13 


33 


<10 








<10 




<10 




<10 




<10 








<1 




13.7 


10.5 


<1 




12.6 


10.8 






<1 




50 


34 


<1 




23 


20 





a In vitro and in vivo estrogen equivalency factors expressed as ngL -1 relative to 17|3-estradiol. 

b References for estrogen equivalency factors: (1) Jobling and Sumpter (1993); (2) Soto et al. (1995); (3) Van Den Belt et al. (2004); (4) Preuss et al. (2006). 

c Estrogen equivalency quotient in ng 17|3-estradiol L _1 calculated by multiplying concentration of compound by the estradiol equivalency factor of the compound and summing. Concentrations reported are the mean±SD 
from 6 to 10 samples collected throughout the experiment. 
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mance. Steroid hormones were isolated by Ci 8 SPE, and eluted 
with 95:5 methanol: water. The residue was dried under N 2 and 
the methoxime:trimethylsilyl derivatives were formed by reaction 
with O-methoxyamine in pyridine followed by bis(trimethylsilyl)- 
trifluoroacetamide with 1 0% trimethylchlorosilane. The derivatives 
were analyzed by GC/MS/MS in splitless mode. The limits-of- 
quantification (LOQ), based on a signal-to-noise ratio of 10, ranged 
from 0.2 to 2ngL _1 . Non-steroid EACs were isolated by CLLE 
from 1-L unfiltered samples using methylene chloride following 
ionic strength and pH adjustment. The extracts were analyzed by 
GC/MS in the full scan and selected ion monitoring (SIM) modes. 
Compound identification was based on matching retention time 
(±0.05 min) and ion ratios (3 ions ±20%) against authentic stan- 
dards. Quantitation was based on an external calibration curve. 
Acidic nonylphenolethoxycarboxylat.es were determined by evap- 
oration to dryness, reaction with 10:90 acetyl chloride/propanol to 
form the propyl esters of the carboxylic acids, and analysis by SIM 
GC/MS (Barber etal., 2000). 

2.6. Statistical analysis and calculations 

Proportional, non-normal scores for degree of nuptial tuber- 
cle prominence, dorsal fat pad prominence, and sperm abundance 
were analyzed by Kruskal-Wallis tests. Remaining data were tested 
for homoscedasticity and analyzed by two-way ANOVA. Concen- 
trations of vitellogenin were log transformed prior to statistical 
analysis. All computations were performed with PRISM 4.0 (Graph- 
Pad Software, Inc., San Diego, CA). Significance was accepted at 
the 5% level. Biological results are expressed as mean ± standard 
error of the mean (SEM). Estrogen equivalency factors (EEFs) and 
estrogen equivalency quotients (EEqs) were calculated as described 
previously (Vajda et al., 2008). 

3. Results 

3.1. Environmental and chemical measurements 

Wastewater discharge from the City of Boulder WWTP dur- 
ing the 2005 and 2006 experiments averaged 0.74 m 3 s _1 (17 
million gallons a day). The mean annual concentrations (Jan- 
uary 2005-December 2006) of ammonia, nitrate, biological oxygen 
demand, dissolved organic carbon, and total suspended solids in the 
WWTP effluent were 6.6, 1 1 .8, 1 4.5, 1 0.0 and 6 mg L _1 , respectively 
(City of Boulder, 2008). 

The City of Boulder WWTP effluent contained a mixture 
of steroidal- and non-steroidal EACs that varied in composi- 
tion, concentration, and potency throughout each experiment 
(Table 1 and Fig. 1). The limits-of-quantifkation (LOQ), based on 
a peak-to-peak signal-to-noise ratio of 10, ranged from 0.2 to 
2ngL _1 for the steroids and 10 to lOOngL -1 for non-steroidal 
EACs. Matrix spikes yielded acceptable recoveries (70-120%), and 
none of the compounds were detected in the blanks. Estrone 
and 17p-estradiol were detected in the WWTP effluent dur- 
ing both experiments (Table 1) at concentrations ranging from 
<2.0 to >75ngL~ 1 . Estriol (1.2±0.7ngL- 1 , n = 7) and the syn- 
thetic estrogen 1 7a-ethynylestradiol (3.4±3.2ngL _1 , n = 7) were 
detected in the WWTP effluent in 2005 (1.2±0.7ngL- 1 , n = 7) 
but not 2006. Steroids were not detected in the REF samples 
in either experiment. Alkylphenolethoxylates and alkylpheno- 
lethoxycarboxylates were detected in all 2005 and 2006 EFF 
samples. Alkylphenolethoxylates concentrations ranged from 450 
to 9800 ngLr 1 and alkylphenolethoxycarboxylates concentrations 
ranged from 890 to 50,000 ngLr 1 ). EFF concentrations of 4- 
nonylphenol declined significantly from 2000 ± 1 300 ng L _1 (n = 7) 
in 2005 to 240±100ngL- 1 (n = 6) in 2006 (p<0.05). When cal- 
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Fig. 1. Estradiol equivalency quotients for wastewater treatment plant (WWTP) 
effluent (EFF) in (A) 2005 and (B) 2006 calculated from measured concentrations of 
individual estrogenic endocrine-active chemicals [17a-ethynylestradiol (EE 2 ), 17(3- 
estradiol (E 2 ), estrone (Ei ) and 4-nonylphenol (4-NP)] multiplied by in vivo fathead 
minnow 17L3-estradiol equivalency factors. Mean estrogen equivalency quotients 
were significantly greater in 2005 than 2006, as were the contributions of 17a- 
ethynylestradiol and 4-nonylphenol (p<0.05). 



culated using in vivo values, there was a significant (p<0.05), 
50% decline in the maximum predicted EFF EEqs from 2005 
(50±34ngL- 1 , n-7) to 2006 (23±20ngL- 1 , n = 6). No such dif- 
ference in EFF EEq was detected when calculated with in vitro EEFs. 

3.2. Biological measurements 

There was no significant effect of exposure or exposure dura- 
tion on survivorship or condition in 2005 or 2006 (p<0.05). The 
response of primary and secondary sex characteristic expression to 
EFF-exposure differed between experiments. EFF-exposure signif- 
icantly suppressed the expression of secondary sex characteristics 
in both experiments. EFF-exposed fish had fewer (p < 0.05) and less 
prominent (Kruskal-Wallis, p< 0.05) nuptial tubercles (Fig. 2A-D) 
and significantly less prominent dorsal fat pads (Kruskal-Wallis, 
p<0.05) (Fig. 2E and F). Primary sex characteristics were signif- 
icantly affected by exposure in 2005, but not in 2006. In 2005, 
GSI was significantly reduced in the 100% EFF-exposed fish on 
day-14 and day-28 (p<0.05). On day-14, reference male GSI was 
significantly greater than initial controls (p<0.05). Histological 
examination of testes revealed that EFF-exposure significantly 
reduced the relative abundance of sperm (Kruskal-Wallis, p < 0.05) 
(Fig. 3A and B). Beginning on day-4 of the 2005 experiment, sperm 
were 'strongly' or Very strongly' expressed in all REF fish at all sam- 
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Fig. 2. (A) Number of nuptial tubercles expressed by male fathead minnows exposed to wastewater treatment plant (WWTP) effluent (100% EFF), Boulder Creek water 
upstream from the WWTP outfall (100% REF), and a mixture of the two during the 2005 male fathead minnow exposure experiment. (B) Number of nuptial tubercles 
expressed by male fathead minnows exposed to 100% EFF, 100% REF, and mixtures of the two during 2006. (C) Prominence of nuptial tubercles expressed by male fathead 
minnows exposed to 100% EFF, 100% REF, and a mixture of the two during 2005. (D) Prominence of nuptial tubercles expressed by male fathead minnows exposed to 100% 
EFF, 100% REF, and mixtures of the two during 2006. (E) Dorsal fat pad prominence expressed by male fathead minnows exposed to 100% EFF, 100% REF, and a mixture of 
the two during 2005. (F) Dorsal fat pad prominence expressed by male fathead minnows exposed to 100% EFF, 100% REF, and mixtures of the two during 2006 [bars with 
different superscript letters or * differ significantly (p < 0.05)]. 



pling dates. Sperm were never Very strongly' expressed in 100% 
EFF or 50:50 EFF: REF exposed fish from day-4 forward. By day- 
14, sperm were 'weakly' expressed or absent in 100% EFF-exposed 
fish. However, expression of testicular sperm abundance was not 
always predictive of secondary sex characteristics. For example, in 
2005 after 28-days of exposure to 100% EFF, the only male fish to 
express both nuptial tubercles and a dorsal fat pad exposure also 
had an absence of sperm. 

Plasma vitellogenin concentrations were significantly elevated 
by EFF-exposure in both experiments (Fig. 4A and B). In the 
2005 experiment, concentration dependence was not observed, 
as vitellogenin concentrations peaked after 7 days among both 
50:50 EFF:REF and 100% EFF exposed males (p<0.05). The 2006 
experiment included more dilute exposures, and concentration 
dependence was observed. In 2006, plasma vitellogenin was sig- 
nificantly elevated by effluent concentrations of 25% and higher 
(p < 0.05). Reference male secondary sex characteristics and plasma 
vitellogenin did not differ from initial controls at any sampling 
period in 2005 or 2006 (p>0.05). HSI was significantly elevated 
by EFF-exposure in 2005 and 2006; on day-28, HSI was signif- 
icantly elevated in the 100% treatment than in REF or 50:50 



EFF:REF treatments. There was no effect of exposure or dura- 
tion on the splenosomatic index in 2005 or 2006 (p<0.05). In 
general, organ indices were less sensitive to effluent exposure 
than secondary sex characteristics, sperm abundance, or plasma 
vitellogenin. 



4. Discussion 

4.1. Environmental measurements and complex mixtures 

The Boulder WWTP effluent contained a complex and dynamic 
mixture of EACs that varied within and between experiments 
(Table 1 ). This estrogenic effluent is discharged to Boulder Creek 
where it accounts for up to 75% of stream flow, which provides habi- 
tat to numerous species and supplies water to downstream users. 
Fish were exposed to dilutions of the WWTP effluent with upstream 
Boulder Creek water under controlled conditions to determine 
whether the reproductive disruption observed in free-living Boul- 
der Creek fish (Vajda et al., 2008; Woodling et al., 2006) could 
be attributed to the WWTP effluent chemistry. Dilutions of EFF 
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Fig. 3. (A) Relative sperm abundance in male fathead minnows exposed to wastewater treatment plant (VWVTP) effluent (100% EFF), Boulder Creek water upstream from the 
WVVTP outfall (100% REF), and a mixture of the two during the 2005 experiment. (B) Histology of REF exposed fish did not differ from initial controls in sperm (S) abundance, 
whereas EFF exposed fish had sperm that was weakly expressed or absent [both photomicrographs are presented at the same magnification]. 



with REF approximated seasonal base-flow conditions and attenu- 
ated downstream concentrations. In Boulder Creek, as with many 
streams in the Arid West, reproduction in most fishes follows a 
prolonged period (August-April) of low in-stream base flows when 
WWTP effluent discharges are most concentrated. 



During the exposure experiments, the WWTP effluent contained 
a complex mixture of steroidal estrogens, non-steroidal estrogens, 
and other wastewater contaminants known to modulate endocrine 
function in vertebrates. Estrone and 1 7p-estradiol were detected in 
the WWTP effluent during both experiments (Table 1 ); estrone and 



1 1 00% Ref ^^ 50:50 Eff:Ref 



100% Eff 







* * % 


it 


* * 




10000. 




















-1 






















E 






















3 
























1000 ^ 




















c 






















§ 






















I 


100! 

10 1 


i 






i 






i 







Day 14 



Day 28 



B 


^H 1 00% Ref V77} 5:95 Eff: Ref I 1 1 0:90 Eff :Ref 


^^25:75Eff:Ref r^^l 50:50 Eff: Ref 1^^100% Eff 




1E7 






c ? c c 
























JL 
















1 000000 1 








b 










E 
1 


100000^ 


ab b 






ab 












c 


10000^. 




rh 


















s 




















































<n 




























o 


10001 


























> 




a a 










i& 




; 












100l 

10- 


1 


1 










1 

















Day 14 



Day 28 



Fig. 4. (A) Plasma vitellogenin concentrations in male fathead minnows exposed to wastewater treatment plant (WWTP) effluent (100% EFF), Boulder Creek water upstream 
from the WWTP outfall (100% REF), and a mixture of the two during the 2005 experiment. (B) Plasma vitellogenin concentrations in male fathead minnows exposed to 
100% EFF, 100% REF, and mixtures of the two during 2006 [REF-exposed males did not differ from initial controls; bars with different superscript letters differ significantly 
(P< 0.05)1. 



AM. Vajda et al. /Aquatic Toxicology 103 (2011)213-221 



219 



the synthetic estrogen 1 7a-ethynylestradiol were detected in the 
WWTP effluent only in 2005, but not 2006. The concentrations of 
steroidal estrogens were consistent with concentrations reported 
for WWTP effluent from trickling filter, or trickling filter combined 
with activated sludge treatment (Johnson et al, 2005). 

Estrogen equivalency quotients (EEqs) (Korner et al., 1999) pro- 
vide a means to consider the chemical mixture effects of the EACs 
(Vajda et al., 2008) and were summed after multiplying the aque- 
ous concentration of 17a-ethynylestradiol, 17p-estradiol, estrone, 
and 4-nonylphenol by the estrogen equivalency factor (EEF) for 
each compound (Jobling and Sumpter, 1993; Preuss et al., 2006; 
Soto et al., 1995; Van Den Belt et al., 2004) using methodology 
described in a previous paper from this group (Vajda et al., 2008). 
EEq values for the WWTP effluent fluctuated over the course of 
each experiment, as did the contributions of individual chemi- 
cals (Table 1 and Fig. 1A and B). Similar short-term changes in 
WWTP effluent composition and estrogen concentration have been 
reported at other sites (Barber et al., 2007; Martinovic et al., 2008), 
and are captured in these on-site, flow-through, in vivo exposure 
experiments. It is important that the experimental design captures 
the variability of EAC concentration because equivalent maximal 
vitellogenin induction can be attained through either intermit- 
tent or sustained exposure to exogenous estrogens (Panter et al., 
2000). 

When calculated using in vivo values, there was a 50% reduction 
in the maximum predicted EEqs in 2006, paralleling observations 
of reduced impacts on fish reproduction, largely due to the absence 
of 17a-ethynylestradiol. The WWTP effluent characteristics were 
complex and the steroidal estrogens were not the only EACs that 
contribute to reproductive effects. In the 2006 experiment, follow- 
ing the publication of the US EPA 4-nonylphenol aquatic water 
quality criteria (U.S. Environmental Protection Agency), there was a 
significant decrease (up to 85%) in concentrations of 4-nonylphenol, 
which also contributes to the reductions in the calculated EEq. 
Similar reductions in 4-nonylphenol concentrations and estro- 
genic effects following regulatory action were observed in the 
United Kingdom (Sheahan et al., 2002). Estimation of EEq inte- 
grates independent, complex analytical chemistry results into a 
biologically relevant index to allow comparisons between loca- 
tions or sampling dates. Although useful in evaluating impacts on 
reproductive endpoints, EEqs have limitations when examining 
complex endpoints such as regulation of behavior where estro- 
gens signal with compound and signaling pathway specific potency 
(Tabb and Blumberg, 2006). This critical complexity is lost when 
the chemical composition and concentration of the WWTP efflu- 
ent are reduced to the single dimension of estrogen equivalency 
quotients. 

4.2. Vitellogen induction, suppression of primary and secondary 
sex characteristics 



(Smith, 1978). Disruption of these secondary sex characteristics 
may have adverse fitness consequences if it alters the outcome of 
territorial contests, decreases access to mates, or alters the qual- 
ity of parental care (Grist et al., 2003; Lange et al., 2001 ; Scott and 
Sloman, 2004). Secondary sex characteristics were demasculinized 
upon exposure to dilutions of wastewater in both experiments, as 
evidenced by a significant reduction in the number and prominence 
of nuptial tubercles (Fig. 2A-D), and in a reduced prominence of 
the dorsal fat pad (Fig. 2E and F). Expression of these secondary sex 
characteristics in male fathead minnows is androgen-dependent 
(Panter et al., 2004; Pawlowski et al., 2004a,b; Smith, 1974) and 
their suppression here is consistent with antagonism by exposure 
to exogenous estrogens (Parrott and Blunt, 2005; Pawlowski et al., 
2004a,b). Although the suppression of secondary sex characteristics 
observed here might reverse upon depuration, even the transient 
suppression of reproductive traits can have adverse fitness conse- 
quences if it coincides with a restricted breeding season. 

Paralleling observations of decreased relative sperm abundance 
in free-living white suckers (Vajda et al., 2008), sperm abundance 
declined rapidly and severely in males exposed to 50:50 EFF:REF 
and 100% EFF. While effects on testicular morphology and sper- 
matogenesis are commonly reported in estrogen-exposed male 
fathead minnows ( Leino et al., 2005 ; Miles-Richardson et al., 1 999a, 
b; Pawlowski et al., 2004a,b), the rapid induction of sperm loss is not 
commonly reported and may be an effect of the complex mixture 
characteristics of the WWTP effluent. Adverse effects of effluent- 
exposure on primary sex characteristic expression was restricted 
to the 2005 experiments, as there was no significant effect on GSI 
or sperm abundance in 2006. 

Induction of vitellogenin is estrogen-dependent (Arukwe and 
Goksoyr, 2003) and elevated levels in male fish are a reliable 
biomarker of exposure to exogenous estrogens (Sumpter and 
Jobling, 1995). Plasma vitellogenin was significantly elevated by 
environmentally relevant mixtures of the WWTP effluent typical 
of in-stream dilutions in both 2005 and 2006 (Fig. 4). Increased 
male plasma vitellogenin concentrations have been correlated with 
increased mortality and decreased fitness in male fathead minnows 
(Thorpe et al., 2007). The elevation in plasma vitellogenin in EFF- 
exposed fish supports the hypothesis that estrogenic contaminants 
are responsible for the reproductive disruption in free-living fish 
downstream from the WWTP outfall (Vajda et al., 2008). 

Reproductive disruption in EFF-exposed male fathead minnows 
was less severe in 2006 than in 2005, but significant effects of expo- 
sure were observed for all biological endpoints with the exception 
that in 2006 there was no significant effect on GSI or sperm abun- 
dance. Although reduced biological affects parallels reductions in 
effluent EAC concentrations, the use of different fish stocks the two 
experiments prevented statistically appropriate comparisons of the 
relative reproductive disruption. 



Evidence of reproductive disruption consistent with exposure 
to EACs was found for each reproductive endpoint examined. 
Male fish began each exposure experiment in stimulated repro- 
ductive condition with mature gonads and prominent secondary 
sex-characteristics, but were rapidly demasculinized after expo- 
sure to EFF and its dilutions. In contrast, male fish in the REF 
treatment remained reproductively stimulated with no signifi- 
cant reduction from initial conditions. The observed reproductive 
responses occurred in the absence of overt toxicity and morbid- 
ity previously reported for fishes downstream of WWTP effluents 
(Tsai, 1973). These results are consistent with exposure to exoge- 
nous estrogenic chemicals (Arukwe and Goksoyr, 2003; Brian et al., 
2005; Kidd et al., 2007; Panter et al., 1 998; Parrott and Blunt, 2005). 

Nuptial tubercles and dorsal fat pads play key roles in male 
fathead minnow competitive, territorial, and spawning behaviors 



4.3. Implications 

The rapid and severe demasculinization of male fathead 
minnows exposed to the WWTP effluent parallels previous obser- 
vations of free-living male white suckers in Boulder Creek (Vajda 
et al., 2008; Woodling et al., 2006) and is consistent with exposure 
to exogenous estrogens. The effects of exposure to dilutions of the 
WWTP effluent did not manifest overt toxicity and mortality, but 
showed rapid disruption of multiple endpoints in fish exposed to 
dilutions that corresponds to the prevalent in-stream effluent con- 
tribution downstream from the WWTP (Barber et al., 2006). This 
integrated application of field and laboratory, biological, chemical 
and hydrological investigative methods to a complex environmen- 
tal system has demonstrated definitive adverse effects of WWTP 
effluent exposure on multiple fish species. 
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